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^j^. Infrared evolution equations for small-x behaviour of the non-singlet struc- 

D . ture functions /(^^ and g^^ are obtained and solved in the next-to-leading 

approximation, to all orders in as, and including running a^ effects. The 
intercepts of these structure functions, i.e. the exponents of the power-like 
small-x behaviour, are calculated. A detailed comparison with the leading 
logarithmic approximation (LLA) and DGLAP is made. We explain why the 
LLA predictions for the small-x dependence of the structure functions may 
be more reliable than the prediction for the Q^ dependence in the range of 
Q2 explored at HERA. 

I. INTRODUCTION 

In recent years the deep inelastic scattering (DIS) of leptons off hadrons has been the 
object of intensive theoretical studies. In QCD the hadronic tensor for electron proton DIS 
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X 



is usually regarded as a convolution of two objects: the "partonic" tensor Wfj,^ describing 
DIS off a (nearly) on-shell parton (quark or gluon) and the probability P^ to find the parton 
in the nucleon. As is well known, they have a different description: whereas the tensor Wfj_iy 
can be calculated within perturbation theory, there are no model-independent methods of 
calculating Pj^, but only the general behaviour with respect to the kinematical variables 
can be predicted. Usually, the hadronic tensor Wf^i, is presented in terms of the structure 
functions (see e.g. [0) /i,2 and gi^2'- 
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where p, m, S are the momentum, mass and polarization of the incoming parton and q is the 
momentum of the virtual photon. The structure functions depend on the Bjorken variable 
X = ~q^ /{2pq), with (1 > x > 0) and Q^ = —q^ > 0. Each of the structure functions 
consists of flavour singlet and flavour non-singlet (NS) contributions. Although the singlet 
contributions accounting for gluon splitting are dominant at small x, the NS contributions 
are also interesting quantities to investigate. They can be expressed in terms of densities of 
quarks (Ag) and antiquarks (Ag) inside a hadron. The non-singlet spin- dependent structure 
function g^^ refers to the combination of (Ag + Ag), whereas the non-singlet spin-averaged 
structure function ff^ is related to (Ag— Ag); combining them one could study the evolution 
of Ag and Ag separately. 

The standard theoretical investigation of DIS structure functions, originally developed 
for the kinematic region of x ~ 1 and large Q^, is made through the DGLAP p| evolution 
equations . This approach accounts for the evolution of the structure functions with respect 
to the photon virtuality Q^, whereas the evolution with respect to x is neglected. In other 
words, the DGLAP accounts for the resummation of all powers of In Q^ and systematically 
neglects the powers of In x, which is obviously correct when x is not small. Despite the fact 
that the DGLAP equations provide a good agreement with experimental data (see e.g. 0), 
they are not expected to work well in the small-x region, where logarithms of x are not 



less important than logarithms of Q^ and therefore must be taken into account to all orders 
in as- This resummation leads to a power-like, or Regge-like, behaviour in x and in Q^ of 
the type f^^ oc {\/Q^ /x)"" . Such a behaviour cannot be obtained within DGLAP, which 



predicts a dependence of the type fi^ oc exp J c\q.{1/x) InlnQ^. 

In perturbative QCD the power-like behaviour was obtained long ago [^ for the singlet 
structure functions of spin-averaged DIS and later in refs. PJ^ for the structure functions 
gi and fi^. However, contrary to the DGLAP, the results of refs. 0-§] are obtained in 
the leading logarithmic approximation (LLA) with a fixed QCD coupling as- Hence the 
exponents calculated in [|-|^ (also called intercepts) for the predicted power-like small- a; 
behaviour of all structure functions contain ag fixed at a scale that is not well defined. The 
very existence of such a scale, sometimes called "a reasonable scale" , and estimates for its 
numerical value have been discussed in many works (see e.g. 0). However, according to the 
Regge theory, the expressions for the intercepts must not contain a^ running with Q^, but, 
on the contrary, they just have to be numbers. 

Thus, in order to obtain realistic values for the intercepts, one has to go beyond the LLA 
and include the QCD running coupling effects in the evolution equations at small x, so that 
the running coupling ag would be integrated out in the solutions. Recently in ref. Q we 
have obtained expressions for the non-singlet structure function^ fi^-, with running QCD 
coupling effects taken into account. The resulting intercept does not explicitly depend on 
as but depends on the infrared cut-off. That calculation prompted us to study the running 
QCD coupling effects for the non-singlet structure functions in more detail. 

The structure function g^^ is similar to fi^., because the same Feynman graphs con- 
tribute to both of them. As a result, they obey the same leading order (LO) DGLAP 
evolution equation. A difference between them arises only in next-to-leading order (NLO) 
corrections to DGLAP and can be extracted from the expressions of the second loop anoma- 
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lous dimensions (see e.g. PJTO[|). 

Concerning the small- a; behaviour, the difference between fj^'^ and g^^ can be related to 
the difference in their signatures: the signature of f{^^ is positive whereas Qi^ has a negative 
one. A difference between them means that some high-order non-ladder graphs contribute 
to fl^^ and Qi^ in a different way. Historically, signatures for scattering amplitudes of 
high energy processes were first introduced in the context of the Regge theory (see e.g. 



Tl[). A forward scattering amplitude M(s, u), where s and u are the Mandelstam variables, 



has positive (negative) signature if it is symmetrical (antisymmetrical) with respect to the 
replacement s <-> m: 



mW(m,s) =mW(s,m), M^'\u,s) = -M^-\s,u). 



(2) 



Basically, any forward scattering amplitude M{s) consists of the two parts: 

M(s, u) = mW + M^'\ M(^) = - (M(s, u) ± M{u, s)) . 



(3) 



In the small-x region, neglecting 0{im? / s) and 0{Q'^/s) terms, one can use u ~ — s. As 
the imaginary part in s, '^sM, for s > 0, corresponds to the cross- sect ion, this cross-section 
may acquire contributions from both M^^^ and M^~^ amplitudes. The partonic tensor W^i, 
of Eq. (|I]) can be considered as the imaginary part of the amplitude M^^, of the forward 
Compton scattering of a virtual photon on a constituent quark: 
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so that 
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and in particular. 
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The amplitude M^y is invariant with respect to the replacement q -^ —q, fx ^^ u which 
corresponds to the crossing transition from the s-channel to the u-channel. However, the 
projection operators multiplying /i and gi behave differently under such a transition: the 
operator ((7^,^ — qfj.qu/q'^) does not change, so the invariant amplitude Mi multiplying it 
does not either. On the contrary, the operator {i^^uXaTnqx/ {pq)) acquires a negative sign, 
so the invariant amplitude M3 must acquire a negative sign too. Therefore, the invariant 
amplitudes 

M(+) = i(Mi(s,n) + Mi(n,s)) (7) 

and 

M(^) = i(M3(s,M)-M3(M,s)) (8) 

have positive and negative signatures, respectively, and the DIS structure functions /i and 
gi can be considered as imaginary parts of scattering amplitudes with positive and negative 
signatures: 



/i = -53,M(+) , g^ = -53,M(-) . (9) 

The first observation of the importance of the contribution of the negative-signature 



amplitudes to the cross-sections was made in the context of QED in ref. |T^. Later on. 



in ref. [|T3[, quark-quark scattering amplitudes with positive and negative signatures were 



calculated in the double-logarithmic approximation (DLA). Using the results of those works, 
it was shown in ref. that in DLA, with as fixed, the intercept of gi^ is larger than the 
intercept of fi^. 

In the present paper we obtain explicit expressions for f^^ and g^^ for small values of x, 
which account for both leading double-logarithmic (DL) and sub-leading single-logarithmic 
(SL) contributions to all orders in QCD coupling, including running Ug effects. As logarithms 
of both X and Q^ are important at small x, we account for both of them, constructing and 
solving a two-dimensional infrared evolution equation (IREE). In order to take running a^ 



effects into account, we use the approach of ref. [T^, with the improvement made in our 
previous work [Q. Also, we discuss similarities and differences between running a^ effects 
in the DGLAP and in our approach. Finally we discuss the region of applicability of the 
LLA and estimate the range of Q^ where it can be safely applied. The paper is organized as 
follows: in Section 2 we derive the small-x evolution equations for the non-singlet structure 
functions and solve them. In Section 3 we obtain the asymptotic behaviour of j^^ and 
g^^ ■ In Section 4 we compare our results with the intercepts of the non-singlet structure 
functions obtained in LLA. In Section 5 our findings are compared with DGLAP predictions 
for the NS anomalous dimensions and we discuss the difference in accounting for running a^ 
effects. Finally, Section 6 contains our conclusions. 

II. SMALL-X EVOLUTION EQUATION FOR THE NON-SINGLET STRUCTURE 

FUNCTIONS 

In the Born approximation, which is the pure NS case, My^y is given by the sum of the 
graphs in Fig. |^. They yield 

M^r = (-9,u + 2x^) 2MS„ + L,^^sA 2Mirl, (10) 

where 



^^■'■Born ~ In 
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with Cq being the electric charge of the incoming quark and s = {p + q)'^, u = {p — qY the 
Mandelstam variables. We have omitted unimportant contributions proportional to g^, qy 
in Eq. (^. 

^Born ^^6 amplitudes of defined signature (cf. Eq. (^). Obviously, only the first term in 
Eq. (|TT]), corresponding to the graph (a) in Fig. [I|, has an imaginary part in s at s > 0, which 
gives the same contribution to the amplitudes of both signatures. According to Eq. (|^) we 
obtain the well-known results: 



jBorn ^ gBorn ^ |^(^ _ ^) _ (^2) 

The property f^°^"' = gBom YQ^aahis true in the next one-loop approximation, with one 
gluon added to graphs in Fig. |[ Only the dressing of graph (a) in Fig. |l| contributes to 
the imaginary part in s at s > 0. This explains why the LO DGLAP splitting functions 
and the LO anomalous dimensions are the same for the structure functions corresponding 
to amplitudes with different signatures. The difference between them arises only at the 
second-loop order. Indeed, let us consider graph (a) of Fig. |^. It is actually a non-ladder 
graph, obtained from the Born graph (b) in Fig. |I|by adding two gluons to it. On the other 
hand, graph (a) in Fig. ^ can be redepicted as graph (b) in Fig. 0. It is easy to check that its 
imaginary part in s (the quark propagators to be cut are marked with crosses in Fig. |^) has, 
in particular, DL terms (~ ln^(l/x)); this graph therefore cannot be neglected. Then, one 
can see that the DL contribution of this graph, which is symmetrical in /i, v (contribution to 
fi^) and the antisymmetrical one (contribution to Qi^) have different signs. This example 
shows that some higher- loop contributions to f^^ and g^^ are different. 

Thus, a possible regular way of calculating f^^ and g^^ consists of the following steps: 

(i) Calculate the forward scattering Compton amplitude M^i, obtained by adding gluon 
propagators to both Born graphs in Fig. |I|. 

(ii) Extract from it the amplitudes Mi and M3 proportional to {—g^u) and ief^uXaQ^p'^ / (pq) 
(cf. Eqs. (^) and (^), which are the positive- and negative- signature amplitudes M^'^^ 
respectively. 

(iii) Calculate (l/7r)S,M(±) . 

So, first we calculate M^,^, and we do it by constructing and solving an IREE for M^j,. 
In order to account for the evolution in both x and Q^ one must consider DL and SL 
contributions to all orders in a^. High order Feynman graphs contributing to M^,^ may 
have both ultraviolet (UV) and infrared (IR) singularities. Whereas UV singularities are 



absorbed by renormalization, IR ones must be regulated explicitly. Providing gluons with 
a mass violates the gauge invariance. In order to save it and to avoid IR singularities, we 
use the Lipatov prescription of compactifying the impact parameter space (see e.g. works 
p!^ , p!3|JT5| ) . In other words, we introduce the infrared cut-off /i in the transverse space (with 



respect to the plane formed by q and p) for integrating over the momenta of virtual particles: 

ki±> II . (13) 

With such a cut-off acting as a mass scale, one can neglect quark masses and still be free from 
IR singularities. As a result the amplitude M^^, depends on the cut-off and we can study its 
evolution with respect to it, thus obtaining the IREE for M^^. Feynman graphs contributing 
to the non-singlet structure functions are obtained from graphs in Fig. |l] by adding gluon 
propagators to them, without breaking the quark lines, however. With logarithmic accuracy, 
the region of integration over transverse momenta of virtual quarks and gluons can be 
regarded for every Feynman graph contributing to M^,^ as a sum of sub-regions so that 
in every sub-region only one virtual particle (quark or gluon) has the minimal transverse 
momentum k±. We call such particles the softest ones although their longitudinal momenta 
can be large (see e.g. []I3|)- K is essential that integrating over k± one includes /i as the 
lowest limit. The softest particle can be either a gluon or a quark. When the softest 
particle is a quark, DL contributions come from the sub-region where another quark has 
the same transverse momentum, so that there appears a t-channel intermediate state with 
a soft quark pair. Therefore in order to keep the DL contributions, one must consider M^,^ 
as convolution of two amplitudes connected by the softest quark pair, each quark with the 
minimal transverse momentum k±, as depicted by graph (b) in Fig. ^ The first amplitude 
appearing in graph (b) (the upper blob) is the Compton forward scattering M^^, and the 
second one (the lower blob) is the forward scattering amplitude of quarks Mq. On the other 
hand, when the particle with the minimal k± is a virtual gluon, it can be factorized out of the 
amplitude with the DL accuracy [fH] . This means, in particular, that the propagator of such 



softest gluon is attached to external lines only, yielding zero result. Having added the Born 
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contribution to the convolution graph (b) in Fig, 3, we arrive at the equation of the Bethe 
- Salpeter type for M^^. As usual the Bethe - Salpeter equation relates M^^ with on-shell 
quarks (Ihs of equation in Fig. ^ to M^j^y with off-shell quarks (rhs of Fig. ^). However, all 
amplitudes in Fig. ^ are actually on-shell because the virtual quark pair in Fig. |^ has the 
minimal transverse momentum fc^. Integrations in blobs imply k^ to be a new IR cut-off 
and, therefore, a new mass scale for these blobs (see ||l3| , p!5| , p!6| ) . 

Let us note that in the DGLAP approach such softest quark pair is always in the lowest 
ladder rung because of the DGLAP ordering: in this approach all ladder transverse momenta 
ki±_ are indeed ordered as 

/i < A;ix < A;2x < ... < Q\ (14) 



the numeration running from the bottom of the ladder to the top. Equation (|l^) shows 
that only the integration over ki±_ has /i as the lower limit. The lower limits for other ki±_ 
with 1^1 are expressed through the other transverse momenta. Therefore the small- a; 
DGLAP equations correspond to the case where the quark scattering amplitude in graph 
(b) in Fig. ^ is considered in the Born approximation. As is well known, the ordering of 
Eq. (|T^ corresponds to accounting for ln"'(Q^//i^) without considering terms ln'"(l/x) not 
accompanied by logarithms of Q^. On the other hand, as we investigate the small-x region, 
we have to lift this ordering, allowing for the transverse momentum of quarks in any ladder 
rung to reach the lowest limit k± = fi. Obviously, this increases the region of integration 
over ki^. Then, a quark pair with minimal k_j_ appears in any rung of the ladder Feynman 
graph. 

The convolution of the amplitudes depicted in graph (b) prompts us to apply the Mellin 
transform for solving this equation. However, to respect the signatures of M^^'> it is more 
convenient to use the asymptotic form of the Sommerfeld - Watson (SW) transform (see 



11 



M(±)( — ,^ 



C^i^ ^-'M^«(4). (-) 



9 



where 

e^^^ = -^-^^ (16) 

is the signature factor, for which this transform differs from that of Melhn. The inverse 
transform to Eq. (|I3]) is 

fWH = -^^/ rfpexp(-a;p)3,MW(p), (17) 

sm(7ra;) Jo 

where p = ln(s//i^). 

The small-x region corresponds effectively to the dominance of the small-ti; region in 
Eq. ([151) . Expanding ^ into series in u and retaining the first terms in Eq. ([15| , ^^'^^ ~ 
— 1 + tnuj/2 and ^'^^^ ^ i7ruj/2, we see that in the small-x (or small-u;) region 



icj,M(-) = i/™|^(i, c.f '*)(.. !^); (18) 



ia.M<-..ir^(4)".f<*.(..«' 
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we need to calculate F^"^^ only to obtain the structure functions ff^^ and g^^ (cf. Eq.(^). 
To do this one can apply the SW transform of Eq. (|1^) to the equation depicted in Fig. [^ and 
then differentiate it with respect to In/i^. On the other hand the Born amplitude does not 
depend on /i and therefore vanishes when differentiated with respect to fi. Then, observing 
that 

~ ^ "dJF' ~ 51n(s/p2) + 51n(QVp2) ( ^) 



corresponds to 



^F(^) + ^ (20) 



for the amplitudes F^'^\ where we have defined y = ln((5^//i^), and differentiating the rhs 
of the equation depicted in Fig. |] with respect to In/i^, we are led to the following IREE : 



^ + u;)F^^\u;,y)=\^^{l + Xu;)]F^^\u;,y)F;i^\u;) , (21) 



w 



where A = 1/2. For more details, see ref. |Q, where the latter equation was obtained for 
F = F^'^\ The only difference now is that IREE (^T]) involves also the negative signature 
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amplitude Fq of forward quark - quark scattering instead of Fq only. The solutions to 



(+) 



Eq. (|TD are 



F(±) = Cexp( 



-LU 



<±) 



Xu)F^^>{uj)/8n' 



y 



(22) 



which contain two unknown quantities. The first one is an arbitrary factor C, essentially of 
a non-perturbative nature, which can be fixed by comparison with the data. 

The other unknown quantity Fq can be specified by constructing and solving the IREE 
for the amplitude Mq of forward scattering of quarks. The equation for Mq is shown in Fig. ^. 
The first two terms in the rhs look similar to the equation in Fig. ^. Indeed, the rhs consists 
of the Born contribution depicted by graph (a), the convolution of two quark scattering 
amplitudes (graph (b)), the intermediate quarks having minimal transverse momenta k±, 
and a new contribution depicted by graphs (c) - (f). For symmetry, the contributions of 
graphs (c) and (d) are equal. The same is true for graphs (e) and (f). All these graphs 
correspond to the case where the particle with the minimal transverse momentum k± is 
a gluon. The propagators of such soft gluons are attached to external lines because by 
integrating over the soft gluon momenta the most singular DL contributions come from the 



region where each soft virtual gluon is factorized out (see e.g. ref. [T^). Self-energy graphs 
are absent in Fig. | because the Feynman gauge is used. The blobs in these graphs imply 
integrations over momenta of internal virtual particles with k± acting as a new IR cut-off. 
Let us consider in detail the contribution of the graphs in the rhs of the equation in Fig. ^. 
The Born contribution coming from graph (a) in Fig. ^ is given by 

5{±) 



-2nCf 
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±as{-s)- 



(23) 



s — /i^ + «e —s — jj,'^ + le 

where the quark mass is substituted by /i in the denominator and is suppressed in the 
numerator as well as —q^. We use for a^ the following formula: 

1 1 



as{s) 



bH-s/klcD) b 
1 



Hs/^QCD 




b[ln\s/AlaD) + ^' ln^(^/A^CD 



7r^ 



(24) 
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where b = (33 — 2n/)/127r (n/ - the number of flavours). 

Equation ( p^ is written in such a way that as{s) has non-zero imaginary part when s 
is positive. Applying the transform Eq. ( p!7D to Eq. (|23| ) we define 

i?(±) = _/ rfpexp(-cup)S,fiW(p) = -, (25) 

vru; JO C(j 

where p = ln{s/fj?), 

7] f°° dpexp{—puj) 






(26) 



r/^ + TT^ Jo (p + r])^ + TT^ 
and we have introduced rj = ln(/i^/AQ^^) assuming p > Aqcd- 

The first term in the square brackets in Eq. (^) corresponds to the imaginary part of 
Eq. ( plD and the second one comes from the imaginary part of l/(s — p^ + te). 
The expression corresponding to graph (b) in Fig. ^ is 

1 



^7r2 



+ A 
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Here, the softest partons (i.e. partons with minimal k±) here are given by the interme- 
diate quark pair. The remaining graphs, (c) - (f ) in rhs of Fig. ^ correspond to the case 
when the softest parton is a gluon and therefore can be factorized out, i.e. its propagator 
can be attached to the external lines in all possible ways. 

The amplitude Mq in the Ihs of Fig. § is colourless, i.e. it belongs to the singlet repre- 
sentation of the group SUc{3). As one of the virtual gluons is removed from the blobs in 
graphs (c) - (f), these blobs are nolonger colourless. They correspond to the colour octet 
scattering amplitude Mg, where k± acts as a new IR cut-off. As Mg does not depend on the 
longitudinal components of k, one can integrate over them, easily arriving at 

3M(t,) = C^ J^^ ^-^<e)^sM^{^y (28) 

Equation (pSj) reads that there is total compensation between the cuts of Mg, so that only 



cuts that do not involve Mg (i.e. two-quark ones) contribute to ^Mf'^_jy 

As we are going to obtain Mg to all orders in a^ with SL accuracy, we must first 
know Mg with the same accuracy. However, for an evaluation of the asymptotic behaviour 
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of Mq , this knowledge is not necessary. Indeed, it was shown in ref. [17| that in DLA 
the blobs in Fig. |^c - f can be approximated, with a few percent accuracy, by their Born 
contribution. The reason is that the octet amplitude Mg decreases with energy so quickly 
that higher-loop contributions can be neglected. Motivated by this result and noticing from 
the results of ref. [§] that, apart from running ag effects, SL contributions do not change 
DL results drastically, we use this approximation in the present work, replacing the blobs 
in graphs (c) - (f ) of the equation in Fig. ^ by their Born value, although with the running 
as- Substituting the Born approximation M^°™ for the colour octet amplitude Mg, 



KMg^ ^ ^M^Born = -2vr3fJ 



aAs] 



s — j2^ + ze 



=Fa5 



— s 



'S — ^^ + le 



-2-K[^as{s)T^as{-s)], 



(29) 



we then obtain 



^ 



M. 



(±) 



^''^^ [^a,{s) - ^a,{-s)] r ^^a^ik'). 



(30) 
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Finally, applying the Mellin transform Eq. ([T7|) to Eq. (|30D , and using Eq. (|2^) for a^, we 
obtain, for contribution of graphs (c) - (f ) in Fig. ^ 
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In DLA, when as is fixed, instead of Eq. (^) one gets 

AalCp 



T 



1 



(p + 77)2 + 7r2 p + ?7 



(31) 
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um 



(32) 



The positive-signature contribution D^~^^ in the DLA is zero, so that Mg does not con- 
tribute to Mq at all. In the Feynman gauge this means that the positive-signature contri- 
bution of non-ladder graphs in the DLA is zero and only ladder graphs must be considered. 



This remarkable observation was first made in ref. ||T2[ for QED processes, and later in |T^ 
for quark scattering. Equation (|3T| ) shows that running ag effects violate the total com- 
pensation of non-ladder graph contributions. In our previous work [^] we had neglected this 
non-compensation. We discuss it in more detail in Section 4. 
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Now we are able to write the full equation for Fq shown in Fig. ^ which generalizes 



the equation for Fq = Fq obtained in our earlier work 



^(±), 



UJ 



A(uj) 1 



1 + A- 



{Ft\u. 
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UJ 



(33) 

UJ »7r^ luj i "• ' ' ' ^ 

where D^^^ and D^~^ are given by Eq. (^); A = 1/2 corresponds to SL contributions not 
related to running coupling effects. 

Equation (^) has the following solution 



F^^ = Att^ 



uj-^juj^-{l + Xu){A{u) + ujD(±){uj))/2n^ 



(34) 



1 + Xuj 

Combining this result with Eqs. (^, (JT^), (^) we finally obtain the following expressions 
for the structure functions f{^^ and Qi^: 



fi 
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(35) 



with C arbitrary. 



III. ASYMPTOTICS OF THE NON-SINGLET STRUCTURE FUNCTIONS 

The small-x asymptotic behaviour of Qi^ and //^"^ can be obtained by evaluating the 
expressions found at the end of last section with the saddle-point method. However, it is 
much easier to get it directly, by noticing that at small x g^^ ~ x^'^o and /{ ~ x~^o , 
where uJq is the rightmost singularity of F^'^\ Equation (p^ reads that this singularity is 
a square root branch point given by a solution toQ 



uj^-{l + Xuj) 



2Cf 



nb J ANCf Jo 



2Cf 
7ib 

oo 



V 



dpe~ 



Ulp 



rf + TT^ Jo (p + riY + 77^ 



dpe-'"' In 



p + T] 



P + V 



+ 



T 
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(36) 



^There is a misprint in the coefficient before the square brackets of Eq. (52) in ref. ^. The true 
coefficient is (2Ci?/7r6). However the numerical results of ref. Q are correct. 
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We recall that p = ln(s/yU^) and r] = \n{p'^/AQQ,j^). This equation can be solved numerically. 
Before doing so, let us notice that if we keep a^ fixed in Eq. (^) and, for self-consistency, 
put A = 0, it transforms into the well-known algebraic equation in the DLA: 



2a.C 
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TT 
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(37) 



with the obvious solutions 
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where the positive-signature leading singularity reads: 
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Now let us come back to Eq. (^). Substituting the value of Fq of Eq. ( ^41) at the point 
UJ = ujq , where the square root turns to zero, we arrive at the power-like asymptotics: 
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(40) 



Obviously, the solutions uJq and cug to Eq. (^) (also called the intercepts of //^ and 
g^^ , respectively) depend on the choice of the parameters n/, /i, Aqcd- We recall that we 
keep fi > Aqcd- Also, fi must be greater than the mass of the heaviest involved quark, 
which fixes nj. In numerical estimates below, we have used Aqcd = 0.1 GeV and nj = 3. 



Equation (pq) has been solved numerically for different values of p. The dependence of the 
intercepts uJq on /i (in terms of 7]) is plotted in Fig. |[ The plot shows that luq' is greater 
than Uq for all fx. Both ujq and uJq first grow with increasing /x, getting maximal at 
fi ^ 1 GeV , where 



OlI, 



(+) 



fi(+) = 0.37 



(41) 



and 
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4"^ = ^^'^ = 0.4, (42) 

and smoothly fall off after that. Equation ( |36D also tells us that contrary to the double- 
logarithmic approximation of Eq. (|38|) , there is no total compensation between the DL 
contributions of the non-ladder graphs to uJq (see expression in the second squared brackets) 
when as is running. That fact was neglected in p]. Once taken into account, it leads to 
a steep increase of both luq and cUq at small rj. However, we regard this increase as an 
artefact of our approach: using expression (^) for a is becoming unreliable at such small rj. 
Another reason for not considering the small-?] region is that when t] < 2.86 the values of 
uJq become complex, which leads to negative f(^^. Thus, we find that our approach is valid 
for 7] > 2.8. Recalling that rj = InKfi/AQCDY], with fi acting as an input for our evolution 
equations, we see that, for self-consistency, we must keep /i > firnin = ^^-^^qcd = 1-75 GeV 
as Aqcd = 0.1 GeV. A further look at Fig. |^ shows that ujq and cOq are maximal at 
1] ^ 5. Then they smoothly (~ I/?]) fall with growing rj. The reason for this fall-off is 
quite clear: at large t], the running ag is approximated by 1/(677). The ?7-dependence of 



LU, 



(±) 
; 



plotted in Fig. |[ has such a complicated form because of the competition between 



the relative weights of the vr^- and //-contributions in the denominators of Eq. (0). In the 
absence of the vr^-terms ujq will depend on r/ in a smoother way, as shown by the curves 
(3) and (4) in Fig. ^j. We will further discuss this plot in the next section. 

IV. COMPARING WITH THE LEADING LOGARITHMIC APPROXIMATION 

The power-like (Regge-like) small-x behaviour of //^'^ and Qi^, Eq. (0), with DL in- 
tercepts ti^o = ujfj^ given by Eqs. (pHD and (|3^), was obtained in refs. and 0. These 
expressions explicitly include the dependence of a^ upon an unknown scale, because the DL 
approximation treats a^ as fixed. On the other hand, the LLA approach postulates that 
asymptotically, at very high energies, LLA results should dominate over all sub-leading con- 
tributions. Having accounted for the SL contributions, we are now able to check the validity 
of the LLA predictions for the NS structure functions. First, let us notice that accounting 
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for running a^ does not lead to intercepts lUq explicitly dependent on a^, although they de- 
pend on Aqcd- Then, thanks to our choice of the infrared regulation, Eq. (|T^) , Uq depend 

-,DL 



on fx. But in spite of such a difference between luq and the DL intercepts cuS.^, one can see 



from Eq. (|3|) and Fig. || that the DL ratio lj^^ - u;^-^ /ujf^J'^ = 0.055 basically holds for 
ujq also, changing slightly from 0.09 at 77 = 3 to 0.06 at 77 = 8. Therefore, the DL result 
for the ratio in Eq. (^) is a good approximation. The vr^-terms that we have accounted for 
in Eq. (PUDare quite important because of their large value. In principle, they are formally 
beyond the SL accuracy we have kept through the paper. Indeed, contributions ~ vr^ may 
appear also from integrations over phase space. Such contributions in higher orders of the 
perturbative series are beyond the control of any kind of logarithmic approximation. On the 
other hand the vr^ we have accounted for appear as the result of respecting the analytical 
properties of the scattering amplitudes. So, we cannot simply neglect them. Moreover, it 
has been shown recently in ref. [0 that a value Vl^^' = 0.4 is in good agreement with the 
experimental data for the structure function F3 . When the vr^ in Eq. ([36| ) are dropped, 
il^~^ = 0.4 corresponds to /i ~ 5.5 GeV (see curves 3 and 4 in Fig. ^. At this value of /i the 
SL contributions not related to the running a^ are small, i.e. one can neglect them, putting 
A = in Eq. (pG]). Therefore we conclude that /iq = 5.5 GeV is a good estimate of the value 
of the mass scale for LLA evolution equations. In other words, yUg acts as a momentum scale 
for the evolution equations obtained in LLA. Therefore these equations are not expected to 
reproduce the correct Q^-dependence until at least 

Q^ > /x^ ^ 30 GeV^ . (43) 

As /ig is pretty close to the typical values of Q^ at HERA, we can see that using the LLA 
evolution equations is not very reliable for such Q^. 

Let us discuss the x-dependence. In refs. [^9]-|21|, it is assumed that a^ depends on Q^ 
in the expression for the intercepts. On the other hand the results of ref. 0] show that there 
is no such dependence. Now we have seen that in LLA a^ should depend rather on /Xg than 
on Q^, and as the Q^ range covered at HERA is not far from /Xq the estimate «« = asiQ"^) 

17 



might still be in reasonable agreement with the data. 



V. COMPARISON WITH DGLAP 



Now let us compare our results with the DGLAP results for the non-singlet structure 
functions. First, let us check that the corresponding DGLAP evolution equations can be 
obtained from Eq. (^). Writing Eq. (pi]) as 
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(44) 



and noting that the term —u in the square brackets is absorbed by changing the Mellin 
factor {s/ii'^Y to x"*^ (see Eq. (|3|)), then 
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(45) 



are the new small-x non-singlet anomalous dimensions. However, when x ~ 1 the main 
contribution of the Mellin factor exp[ct; ln(l/a;)] comes from the region of rather large values 
of tu, where u;ln(l/x) < 1. That makes it possible to expand the exponent Fq in Eq. (pSj) 
into a series in l/uo. Doing so, we obtain 
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(46) 



Retaining only the first term in the rhs of Eq. (^61), one arrives at the expression 



~(±)_ ^ 



1 + A 



(47) 



for the leading order anomalous dimension, which is similar to the expression for the leading 
order DGLAP non-singlet anomalous dimension 
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(48) 



Let us first compare the two results by assuming that a^ is fixed. Equations (H^ and 
f) coincide completely because in this case A = AirasCp. Let us then compare the next 



two terms in Eq. (|^) with the second order DGLAP anomalous dimension '^(2) , again for 
fixed as- Retaining the leading terms, proportional to 1/u^ , we can see that 
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transform into (see Eq. (p2| 
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cu3 



which coincide with the leading terms of Eq. (B18) of ref. |]I0|, also used in ref. for 
the second order non-singlet anomalous dimensions. Unfortunately, we do not obtain exact 
coincidence between non- leading terms, which are proportional to l/cu'^, but this can be 
explained by a different choice of the factorization and regularization procedures. Therefore, 
the non-singlet anomalous dimensions given by Eq. (^) correspond to a resummation of 
7(1) to all orders in the QCD coupling and to accounting for the signature-dependent DL 
contributions to order a^. Such a resummation leads to the power-like small-x behaviour, 
which cannot be achieved by incorporating any finite number of NLO contributions into 
expressions for DGLAP anomalous dimensions. 

Another difference between TREE and DGLAP is the different treatment of the QCD 
coupling. Although both approaches use the same formula (Eq. (p4D) for a^, the DGLAP 
uses as = as{kfj^ in i-th ladder rung, ki± being the ladder quark transverse momentum 
whereas TREE suggests that in the i-th rung a should depend on the gluon virtuality [| 
{ki — ki-iY ■ Let us explain this difference. When a; 7^ 1, the Born contribution to f^^ is 
zero and f^^ obeys the Beth - -Salpeter equation 



^We use here the standard numeration for the momenta of the virtual ladder quarks fcj, so that 
i = corresponds to A;o = p and, increasing the number of runs, it leads to the top of the ladder. 
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f^s = fd^kmq + k) 



2^ 1 ^ asijp-ky) 



(52) 



We have suppressed all unimportant factors in Eq. (|5^) and have used the notation k instead 
of ki for the ladder quark momentum in the lowest rung. We have noted as $ for off-shell 
f^^ . Usually it is convenient to integrate over k in (|5^), using the Sudakov parametrization: 

(53) 



k„ = a{q + xP)y + (3py + fcj 



so that 



2{pk) = 2{pq)a, 2{qk) = 2{pq){(3 - xa), k'^ = 2{pq)ap - k] 



But here it is more convenient to use the variable 



(54) 



m^ = (p — ky 



(55) 



instead of the Sudakov variable a, writing Eq. (|5^) as 
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The argument of $ in Eqs. ( |5^ ) and (^) should be positive, i.e. 
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When a; ~ 1, one can neglect the last term in Eq. (|57|) . It leads to (/5 — a;)(s(l — /3) —m^) > 
which allows m^ to be neglected. Therefore, instead of Eq. (|56D , we obtain 



f^s ^ f d^k^df3dm^<!> {[s{(3 - x){l - (3)] (1 - (3)-^) 53- 



(58) 



[/?m2 + A;2_][m2] ' 

After that, one can apply the Cauchy theorem relating an analytic function to its imaginary 
part. This permits to perform an integration over m^ by taking residues. As m^ is positive, 
the only residue is m"^ = —k'^/j3. Thus, we obtain 
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(59) 
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As /3 > X and x ~ 1, one can drop f3 in the argument of as and arrive at 



which corresponds to the DGLAP. Usually the negative sign of the argument of a^ is dropped 
together with the negative sign of the argument in Eq. (p^. As the argument of a^ in Eq. ( |59D 
is negative, the expression for a^ in the DGLAP approximation of Eq. (|60|) does not contain 
vr^-terms at all, while the argument of as in Eq. (|5^) is positive and therefore there are 
vr^-terms in Eq. (^2]). 



Obviously, the transition from Eq. (^) to Eq. (|59| ) holds only when a; ~ 1, which is 
the DGLAP kinematics, and fails for small x. In particular, neglecting kj_ in Eq. (|57D 
has an important consequence for the DGLAP evolution, where one picks up logarithmic 
contributions through integrations over k±. At small x one should use the ttg-dependence 
given by Eq. (^) rather than that given by Eq. (|60|). 

VI. CONCLUSIONS 

In conclusion, we have obtained a generalization of the second-order DGLAP evolution 
equations for the non-singlet structure functions at small x. Besides resumming the LO 
DGLAP anomalous dimension to all orders in ag and accounting for running as effects, we 
account for the difference between ff^ and Qi^, which is due to the difference in the second- 
order anomalous dimensions. We have shown that, with single- logarithmic contributions and 
running as taken into account, ff^ and g^^ have the scaling-like asymptotic behaviour 

as also obtained earlier in the DLA [^,|| at asymptotically small x. We have appropriately 
calculated the intercepts uJq . We have made detailed comparisons of our results with the 
predictions obtained in the leading-logarithmic approximation and with the DGLAP. Results 
obtained with the evolution equations usually depend on a mass scale nl acting as a starting 
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point of the evolution with respect to Q^, and we estimate the value for this mass scale for a 
LLA evolution equation as /Iq ~ 30 GeV^. This implies that these equations can reproduce 
the correct Q^-dependence for the non-singlet structure functions (and likely, for the singlet 
ones) only for approximately Q^ > 30 GeV^ . On the other hand, as far as the x-dependence 
is concerned, we have also shown that the estimate a^ = as(Q^) used in refs. |]I^|-pT|] in LLA 
can be in a reasonable agreement with the x-behaviour of the structure functions for the 
HERA data. 

VII. ACKNOWLEDGEMENTS 

The work is supported in part by grant INTAS-97-30494, and by EU QCDNET contract 
FMRX-CT98-0194. 



22 



REFERENCES 

[1] R.P.Feynman, Photon-hadron interactions (Benjamin, 1972); B.L.Ioffe, V.A.Khoze and 
L.N.Lipatov, Hard Processes (North Holland, 1984). 

[2] G.Altarelli and G.Parisi, Nucl. Pliys. B126 (1977) 297; V.N.Gribov and L.N.Lipatov, 
Sov. J. Nucl. Phys. 15 (1972) 438; L.N.Lipatov, Sov. J. Nucl. Pliys. 20 (1972) 95; 
Yu.L.Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[3] G.AltareUi, R.D.Ball, S.Forte and G.Ridolfi, Nucl. Phys. B496 (1997) 337. 

[4] V.S.Fadin, E.A.Kuraev and L.N.Lipatov, Sov. Phys. JETP 44 (1976) 443 and 45 (1977) 
199; Y.Y.Balitskij and L.N.Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

[5] B.I.Ermolaev, S.I.Manayenkov and M.G.Ryskin, Z. Phys. C69 (1996) 259. 

[6] J.Bartels, B.I.Ermolaev and M.G.Ryskin, Z. Phys. C70 (1996) 273 and C72 (1996) 627. 

[7] S.J.Brodsky, V.S.Fadin, V.T.Kim, L.N.Lipatov and G.B.Pivovarov, JETP Lett. 70 
(1999) 155-160; |hep-ph/990l22g. 



B.I.Ermolaev, M.Greco and S.I.Troyan, Nucl. Phys. B571 (2000) 137; |hep-ph/9906'77B 



[9] M.Gluck, E.Reya, M.Stratmann and W.Vogelsang, Phys.Rev. D53 (1996) 4775. 

[10] E.G.Floratos, C.Kounnas and R. Lacaze, Nucl. Phys. B 192 (1981) 417. 

[11] P.D. B.Collins, An introduction to Regge theory and high energy physics. Cam- 
bridge, 1977. 

[12] V.G.Gorshkov, L.N.Lipatov and M.M.Nesterov, Yad. Fiz. 9 (1969) 1321. 

[13] R.Kirschner and L.N.Lipatov, Nucl. Phys. B213 (1983) 122. 

[14] L.N.Lipatov, Phys. Lett. B116 (1982) 411. 

[15] B.I.Ermolaev and L.N.Lipatov, Int. J. Mod. Phys. A4 (1989) 3147. 



23 



[16] B.I.Ermolaev, V.S.Fadin and L.N.Lipatov, Sov. J .Phys. 45 (1987) 508. 

[17] B.I.Ermolaev and S.I.Troyan, in Proceedings of 5th Intern. Workshop on DIS, New York 
1998, p. 861. 

[18] A.Kataev, G.Parente and A.V.Sidorov. Nucl. Phys. A666/667 (2000) 184; Nucl. Phys. 
B573 (2000) 405. 

[19] B.BadaIek and J.Kwiecinski, Phys. Lett. B418 (1998) 229. 



[20] J.Kwiecinski and B.Ziaja, |hep-ph/9902440 . 



[21] Y.Kiyo and J.Kodaira, |hep-ph/ 98034481 |hep-ph/9711260| ; Z. Phys. C74 (1997) 631 



24 



FIGURES 
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FIG. 1. The Born graphs for the DIS ampUtude M^,y. 
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FIG. 2. Two-foop graph obtained from Fig. ^, contributing to QsM^jy (a), and the corre- 
sponding physical process (b). 
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FIG. 3. The evolution equation for the DIS amplitude M, 



/lU' 



25 
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FIG. 4. The evolution equation for the quark scattering amplitude Mq. 
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FIG. 5. Dependence of the intercept loq on infrared cut-off r] = ln{fj? /Aqcd)- 1: for fi , 2: 
for gi'^, 3 and 4: for f^^^ and gi'^, respectively, without accounting for vr^-terms. 
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